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Abstract 



The total and differential cross sections for the production of triply charmed Occc baryons in e~^e 
annihilation are calculated at the Z-boson pole. 
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Oi 1- Introduction 

Q ■ Investigation into the properties of baryons containing two or three heavy c and b quarks, 

i-C ' the features of their production at operating accelerators and those under construction, and 

ji^! their lifetimes and decay modes is topical in particle physics, but these issues have not 

^ ■ yet received adequate study. All that is currently known in these realms from experiments 

amounts to the claim [1] that a doubly charmed baryon was observed in experiments with 

^ . a beam of charged hyperons at FERMILAB. Theoretical investigations of baryons containing 

d • two heavy quarks are reviewed, for example, in [2]. Calculations available in the literature 

that deal with the cross sections for the production of baryons containing two heavy quarks 

treat primarily processes described in the fourth order of standard perturbation theory-that 

is, processes leading to the production of respective diquarks [3]. Only in [4] were sixth-order 

calculations performed, where the process e^e^ —>■ ssccbb was associated with the production 

of an Qgcb baryon in e~^e~ collisions. The production of baryons involving three heavy quarks 

has not yet been considered. 

The present article reports on a continuation of the investigation begun in [4] , providing 

a description of some features of the process involving the production of triply charmed 

baryons Qccc in e^e~ annihilation. This case has nothing to do with the production cc 
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diquarks, because they can transform, with a probabiUty close to unity, only into 5++ or 
"ccd baryons, a negligible fraction of these diquarks going over to Qccc baryons. As a matter 
of fact, calculations in the sixth order of perturbation theory for the elementary process 
e+e~ — > cccccc are the only possibility of theoretically studying triply charmed baryons. 
The main contribution to the amplitude of this process comes from 504 Feynman diagrams. 
In relation to the production of Q,scb baryons, which was considered previously and where all 
components have different flavors, the calculations for Qccc baryons are complicated by the 
need for taking into account the interference between identical particles. 

In studying the production of Qccc baryons in proton-proton collisions, it would be neces- 
sary to consider an order of magnitude greater number of Feynman diagrams corresponding 
to the subprocesses qq — > cccccc and gg — > cccccc. Moreover, the description of baryon pro- 
duction in hadron-hadron collisions would require a much greater effort in calculating the 
contributions to the amplitude of the production process from various parton color states 
than in the case of e'^e~ annihilation. 

In present study, the uniflcation of three charmed quarks into an flccc baryon is de- 
scribed within the well-known nonrelativistic approximation [5]. Upon obtaining numerical 
results for the cross sections describing Q^cc baryon production, we analyze the possibility 
of constructing their approximate analytic description in terms of one known fragmentation 
function or another. 



2. Amplitude of flccc production in e+e annihilation 

We assume that the amplitude of the production of triply charmed baryons flccc in e+e" 
annihilation corresponds to the elementary process 

e+(A;i) + e~{k2) c(pi, + c(p2, ^2) + c(p3, ^2) + c(p4, Xi) + c(p5, X2) + c{pe, xs), (1) 

where ki and k2 are the 4-momenta of colliding particles; pi, ...,Pq are the 4-momenta of 
product partons; and and Xj (^^ J = 1; 2, 3) are the color indices of quarks and antiquarks, 
respectively. As usual, we disregard the contribution of the electroweak interaction of quarks 
to the amplitude of process (1), since it is an order of magnitude less than the corresponding 
contribution of QCD interaction. Thereupon, all Feynman diagrams to be taken into account 
for process (1) reduce to the nine basic diagrams in Fig.l, which correspond to different 
positions of the quark-gluon vertices. Thirty-six nonequivalent dispositions of quark and 
antiquark lines characterized by speciflc 4-momenta, polarizations, and color indices are 
possible for each of the basic diagrams 1-7, and 18 nonequivalent dispositions of such lines 
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are possible for the basic diagrams 8 and 9. Since a collision between an electron and 
a positron leads to annihilation into either a photon or a Z boson, the total number of 
relevant Feynman diagrams is 576. 

First, we consider the color structure of the amplitude of Qccc baryon production. Since 
an electron, a positron, and any baryon are singlets with respect to the SU{3)c color group, 
three product antiquaries c must also form an 5'C/(3)c-singlet state. Therefore, the final state 
of process (1) must be fully antisymmetric in the color indices of the three charmed quarks 
bound into an ficcc baryon and in the color indices of the three product charmed antiquarks. 
This requirement, together with the requirement of an appropriate normalization, is satisfied 
by introducing, in the amplitude of process (1), the product (e^'^^^^^ / \^){e^'^^'^^^ /\/6) of 
antisymmetric tensors and performing summation over the color indices of and Xj {hj — 
1,2,3). 

We set = A"/2, where A" (a = 1, ...,8) are the Gell-Mann matrices, and denote by 
N the total number of such permutations of different pairs of color indices of quarks and 
antiquarks that transform the sets (^1,^2,^3) ^-iid (xi,X2,X3) into the sets (61,62)63) ^-i^d 
iXhiXi'iiXki)^ respectively. The color factors associated with diagrams of the types -Z-7can 
then be found by means of direct analytic calculations. The result is 

a,b,C Ci,6.53 Xi:X2,X3 " 

In the sum in expression (2), the index 6) appears twice (as it must)-directly in the 
tensor £?i?2& ^j^j indirectly as the substitute of one of the indices 61 ) 62 ^^^^ 63- The same 
is true for the other Greek indices in the above sum, with the exception of and for the 
zeroth color factors corresponding to diagrams of types 8 and 5, for which we have 

^ ^ ^ •' SilXji ?i2Xj2 ?»3XJ3 ^ ' 

a,b,c 6,6,6 Xi,X2,X3 

where /"^'^ are the structure constants of the SU{3) group. The proof of the equality in 
(3) is given in [4]. This equality means that the total contribution to the amplitude of the 
process in (1) from the diagrams involving three-gluon vertices vanishes. Thus, the number 
of contributing diagrams reduces to 504. 

Since the contribution to the amplitude of process (1) from the diagram that differs from 
a specific diagram by a permutation of N fermion pairs involves the Feynman factor (—1)^, 
it can be concluded, with allowance for (2), that all terms appearing in the amplitude of ftccc 
baryon production have the same sign. 

It should be noted that, in our calculations, we used an additional simplifying approxi- 
mation, setting the c-quark mass to zero in all expressions entering into the numerators of 
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fermion propagators and in all traces. At the same time, we set rric — 1.5 GeV and pj = 2.25 
GeV^ (i = 1,...,6) in all of the denominators of the propagators of virtual particles and 
in the expression for the final-state phase space (of course, the amplitude containing zero 
c-quark mass in the denominator would diverge). But if we used a nonzero c-quark mass 
everywhere in the amplitude and in the square of the relevant matrix element, the volume of 
information to be saved in the computer memory and the time required for the compilation 
of codes and for numerical calculations of the cross sections would grow enormously, which 
would render the problem in question unsolvable with our means. 

In order to estimate the effect of the above approximation on the accuracy of the nu- 
merical results, we repeated the calculation of the cross section for Qsc& baryon production 
in a similar approximation and compared the results obtained in this way with the results 
of the full calculation performed in [4]. It turned out that the cross sections obtained for 
D,scb baryon production within the " massive" and " massless" (for all quarks simultaneously) 
approximations differ only by 8%. It seems reasonable to expect an inaccuracy on the same 
order of magnitude for flccc baryon production as well. Anyway, this inaccuracy does not 
exceed other theoretical uncertainties associated, for example, with the choice of the renor- 
malization scale in the strong-interaction coupling constant or with the wave function for the 
triply charmed heavy baryon. Thus, this approximation appears to be numerically justified. 

Taking into account the aforesaid, we can represent the matrix element for process (1) 
in the form 

D'-^D'' (4) 

Qcos^ 0w{s - M^ + iMzTz) 9s ' ^ ^ 

where 

^ y: e {[{pj+Pi+Pi'r-mmh+k2-Pk'r-mi]-'x (5) 

i,j,fcS{l,2,3} i',j',fc'e{l,2,3} 

X {pi + Pi')~^{pi + Pj + Pi' + Pji)~'^u{pj)Y{Pj +Pi+ Pi'hsv{-Pf) X 

X u{Pkh\h + k2- Pk'hei9v - 9Al5)vi-Pk') + 

+ [iPf +Pi+ Pi'f - mlY^h + k2- Pkf - ml]-^{pi + Pi/)-2 X 

X {Pi+Pj +Pi' + Pj')~'^u{pj)-fs{-pj' -Pi-Pi')Yt^{-Pj') X 

X u(pk)^e(-ki - k +Pkh\9v - 9A'y5)v(-pk') + 

+ [{Pj' + Pt + Pt'Y - ml]-'^[{ki + k2-pk'f - ml]-\pi + pi')-"^ x 

X {Pi + Pj + Pi' + Pj')~^u{pj)-f5{-Pj' - pi- Pi'h^vi-Pj') X 

X u{pk)-f\ki + k2- Pk')le{9v - 9Alb)'v{-Pk') + 
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+ [{Pj+Pi+Pi')^ - ml]-^[{ki + k2-Pkf - ml]-^{pi+Pi,)-'^ x 

X iPi + Pj + Pi' + Pj')~'^^{V>j)l''{P3 + Pi + Pi>)l5v{-\>j') X 
X U{^k)le{-kl - h +Pkh\9v - 9Al5)vi-Pk') + 

+ [{Pj + Pj'+Pkf - ml]-'^[{ki + k2-pk'f - m^^~^{pi + pj/)"^ x 
X {Pj + Pj')~^'"(Pi)75^(-Pi')^(Pfe)7'^(Pj + Pj' + Pfc) X 

X Yih +k2- Pk'heiQv - 9Al5)v{-Pk') + 

+ [{Pj +Pj' +Pkf - ml]~^[{pi+pi> +Pk'f - ml]~^{pi+pi')~'^ X 

X (Pj + Pj')~'^u{Pjh5v{-Pf)u{pk)-f^{pj + Pj' +Pk) X 
X le{-Pi - Pi' - Pk')Y{9v - 9A15)V{-Pk') + 

+ [{Pi + Pi' +Pk'f - m^r^[(ki + k2 - Pkf - ml]-~\Pi+Pi')~^ X 
X {Pj +Pj')~'^u{pj)-fsv{-pj>)u{pk)%{-ki - h+Pk) X 

X iH-Pi - Pi' - Pk')Y{9v - 9'kl5)v{-Pk')} X 
X u{pi)juv{-pi')v{-ki)j%glr - c/^75)M(k2), 

while the expression for can be derived from by means of the substitutions giy — > 1, 
9a 0) ffy ^ Qc = 2/3 and 0. Summation in (5) corresponds to 36 permutations of 

quark and antiquark hnes in diagrams of types 1-7. 

3. Method of orthogonal amplitudes 

In order to derive the expression that is obtained for the square of the matrix element 
upon summation over the final-fcrmion polarizations and averaging over the polariza- 
tions of colliding particles, we use the method of orthogonal amplitudes and the REDUCE 
computer system for analytic calculations. The method of orthogonal amplitudes was pro- 
posed in [6] and was employed in calculations referring to flscb baryon production in e'^e~ 
collisions [4]. 

A simple and mathematically rigorous validation of the method of orthogonal amplitudes 
is the following (to the best of our knowledge, it has not yet been given anywhere). Suppose 
that four-component spinors u{p) and u{p') describing particles of mass m and m', respec- 
tively, their 4- momenta being p and p' {p^ — m^, p'^ — m'^), obey the Dirac equation. Of 
the four linear homogeneous equations for the components of the spinor u{p) [u{p')], only 
two are independent; therefore, each of the four components under consideration can be 
represented as a linear combination of two arbitrary independent constants, denoted here 
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by X and Y {X' and Y'). Any quantity of the form ?7(p')i?-u(p), where R is an opera- 
tor specified in terms of the 7-matrices and their contractions with some 4-vcctors, can be 
represented as a hnear combination of four independent elements XX'*^ XY'*, YX'*, and 
YY'*. Therefore, quantities of the form u{p')Ru{p) can be treated as vectors of a hnear 
four-dimensional space L spanned by the above elements. Any four linearly independent 
quantities of the form Wn = u{p')Onu{p), where the operator On is either unity, 75, V, 
or {V"V"' - V"'V")/2 (with V, V, V", and V" being arbitrary 4- vectors), can be taken for 
basis vectors of the space L. The scalar product {wn, Wn') of vectors Wn and Wn' belonging to 
the linear space L is defined as the product WnW^, summed over the polarizations of fermions 
that are described by the spinors u{p) and u{p'). 

Here, we take, for basis vectors of the space L, four quantities Wn specified by the op- 
erators Oi — 1, O2 — K, Os — Q, and O4 — KQ with the 4- vectors K and Q here being 
orthogonal to the 4-momenta p and p' and to each other - that is, K^p^ — 0, Kfj_p'^^ — 0, 
Q/ip'' = 0, Q^^p'^ = 0, and K^Q^ = 0. Otherwise, the 4-vectors K and Q are arbitrary. They 
can be specified, for example, by the relations K^^ — e^^^'^p^p'^Ua and ~ e'^^'^pyp'pKfj, 
where the 4-vector 0^. is entirely arbitrary. From the orthogonality of the 4-vectors K and 
Q, it follows that KQ = {KQ — QK)/2. The four quantities Wn used are orthogonal to 
one another, {wn,Wn') = CnSnn', Cn 7^ 0, this proving their linear independence and justi- 
fying the name "orthogonal amplitudes." Thus, it was shown that any quantity of the form 
'u(p')-R-u(p) can be represented as a linear combination of orthogonal amplitudes. 

In order to solve the problem of calculating the square of the matrix element, we first 
introduce basic orthogonal amplitudes as 

Wn = U{pi)v{-p3+i), Wi2 u{pi)kiv{-p3+i), (6) 

Wi3 = u{pi)QiV{-p3+i), Wn = u{pi)kiQiv{-p3+i), 
Wei = t^(-ki)M(k2), = v{-'ki)keU{]<i2) , 

Wes = v{-ki)Qeu{k2), w^^ = v{-'ki)k^Q^u{\i2) , 

where i = 1,2,3. We would like to note that the pair combinations of the spinors u{pi) and 
v{—pj) can be chosen in six equivalent ways. 

On the basis of the quantities in (6), we construct 256 orthogonal amplitudes as 

Wnrst = WinW2r'W3sWet, (7) 

where n, r,s,t — 1, 2, 3, 4. 
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The expansion of the matrix element (4) in the amphtudes given by (7) has the form 

4 

= (^nrstWnrst- (8) 

n,r,s,t=l 

In order to derive the coefficients Cnrst in this expansion, we multiply both sides of (8) 
by the factor Wn'r's't', take the sum of the result over the polarizations of all of the fermions, 
and make use of the orthogonality of different amplitudes. As a result, we arrive at 

Cnrst = { X] •^'tJ^nrst}/{'^nrst,Wnrst), (9) 
polar. 

where (wnrst^Wnrst) IS an analog of the scalar product defined above in the linear space L - 
that is, the sum of the squared modulus of the amplitude Wnrst over the polarization of all 
fermions. Since Wnrst involves arbitrariness associated with the choice of the 4- vectors K and 
Q in (6), there is also arbitrariness in the coefficients Cnrst in (9). The substitution of these 
coefficients into (8) leads to an identity whose left-hand side is determined unambiguously. 
Thus, summation on the right-hand side of (8) removes the above ambiguity. 

Since electrons and positrons are treated as massless particles and since the charmed- 
quark mass is set to zero in the numerators of each term of the amplitude for process (1) 
and in respective traces, it is clear that 192 of the 256 coefficients in expansion (8) vanish, 
because they are linear combinations of the traces of an odd number of the Dirac 7 matrices. 
We further hst 64 orthogonal amphtudes in formula (7) that generate nonzero expansion 
coefficients: t = 2,3 with either n,r,s = 2, 3, or one of the indices n, r, s is equal to 2 or 3, 
while the other two belong to the set {1,4}. 

It can clearly be seen that the expression obtained for the square of the matrix element 
upon summation over the polarizations of final fermions and averaging over the polarizations 
of initial particles takes the form 

\M\^^'^ \Cnrst\'^ ■ {Wnrst, Wnrst)- (10) 

n,r,s,t 

In actual calculations by the method of orthogonal amphtudes, we compose one RE- 
DUCE code for traces and tensor contractions that corresponds to 504 terms in any quantity 
M.w^j.gf and then, by means of any text editor (for example, "joe"), perform obvious changes 
necessary for obtaining the REDUCE code for calculating all 64 nonzero quantities A4w^j.g^. 



4. Cross sections for Qccc baryon production at the Z 
pole in e+e~ annihilation 
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In describing the ficcc baryon as a bound state of three charmed quarks, we use the 
nonrelativistic approximation [5]. This means that we disregard the relative velocities of the 
c quarks confined within the baryon - that is, in the laboratory frame, the velocities and 
momenta of all three c quarks produced in process (1) are taken to be identical and equal to 
one-third of the momentum p of the ftccc baryon having the mass M—Snic- With allowance 
for the unification of three charmed quarks into the baryon, the phase space of process (1) 
effectively becomes the 4-particle phase space of the process 

e+(A;i) + e-{k2) ^ Jlccc(p) + c{P^) + c{P5) + c{Pe)- (H) 
The differential cross section for process (11) takes the form 



(27r)^|M|2 |^/;(0)| 



2 



d(7= TT^S {ki + k2-p-p4-p5-p6) X (12) 



X- 



2s M2 
d^p d^pi d^p5 d^pe 



(27r)32£; (27r)32£;4 {27r)^2E5 {2t:Y2Eq 
where ■0(0) is the value that the respective wave function takes in the case where all three c 
quarks forming the ficcc baryon are located at the same point, so that their relative coordi- 
nates are zero. The numerical value of |'0(O)p is taken to be identical to that in [7], where 
it was 

|0(O)p = 0.36-10-^ GeV^ (13) 

In calculating the total and differential cross sections, wc employed codes for numerical 
integration that are based on the Monte Carlo method and which are contained in the 
CompHEP package [8], which is broader. It appeared that the maximum computational 
errors in the differential cross sections came from the first iteration. Therefore, only the 
total cross section was calculated in the first iteration, while both the total cross section 
and the differential cross sections were determined in the next five iterations. Each iteration 
involved 200 000 Monte Carlo calls on the integrand. The errors in calculating the total 
cross section amounted to 1.0%, while the errors in calculating the differential cross sections 
were predominantly 2 to 3% (this is reflected below in the text and in the figures). As was 
indicated above, the error associated with the disregard of the charmed-quark mass in the 
numerators of the amplitude for process (1) and in the traces is a few percent. Moreover, 
we additionally tested the consistency of the cross-section values for two different admissible 
choices of the 4- vectors K^, Q^, Kf, and Qf (i = 1, 2, 3) used to construct the quantities in 
(6), which specify the orthogonal amplitudes (7). 
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In addition to statistical errors, the calculations contain unavoidable theoretical uncer- 
tainties. First, there is the uncertainty associated with the running strong-interaction cou- 
pling constant as a function of the renormalization scale. Since all of the calculations were 
performed at an energy value that corresponds to the Z-boson pole {^/s — 91.2 GeV), it is 
reasonable to specify the couphng-constant values as follows: — as{Mz/2) — 0.134 and a 
= a{Mz) = 1/128.0; accordingly, sin^6'vK = sin^OwiMz) = 0.2240. However, it is not evident 
why it is Mz, and not, for example, the invariant mass of some product quark pair or even 
the Qccc baryon mass, that should be chosen for the characteristic scale of strong interaction. 
Since the cross section for process (11) is proportional to the fourth power of the strong- 
interaction coupling constant, this source of errors is the most important. Second, the ac- 
curacy of the potential model employed as a basis for calculating the baryon-wave-function 
value -0(0) is uncertain. 

For the chosen set of model parameters, the total cross section for the process (Jtot and 
the forward (backward) production asymmetry at the Z-boson pole are 



where apicrB) is the cross section for the production of an Qccc baryon moving in the forward 
(backward) direction with respect to the electron-momentum direction. The cross-section 
value in (14) is close to that of the total cross section for Q,scb baryon production in e'^e~ 
collisions (0.0534 ± 0.0014 fb) if the strange-quark mass is set to 300 MeV [4]. 

The differential cross sections with respect to the transverse momentum px and the 
rapidity Y of Qccc baryons are presented in Fig. 2. The distribution da/dY peaks at a 
small positive value of Y, while da/dpr has a maximum at pr ~ 12 GeV. Note that the 
maximum of the differential cross section with respect to the transverse momentum of Qccc 
baryons occurs at a Pt value much lower than that for Qgcb baryons produced under the 
same conditions, in which case da/dpx peaks within the Pt interval 23-26 GeV. 

It is desirable to associate our numerical results with some simple analytic form, which we 
will seek among well-known fragmentation functions [9-12]. It is clear that the production of 
a triply charmed baryon can hardly be interpreted as a fragmentation process, since each of 
the three c quarks can be treated, on equal footing, as a fragmenting quark produced at the 
^/Z vertex and since the interference between identical quarks is likely to be significant in the 
process being considered. However, we accept not the physical concept of the fragmentation 
model but its mathematical form used in processing experimental data on e+e" annihilation 



(7tot = (0.0404 ± 0.0004) fb. 



(14) 



Afb = {ctf - crB)/{crF + (Jb) = 0.101 ± 0.005, 



(15) 
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(see, for example, [13]); namely, we set 

^ = (Tec • L>c^n,,,(^), (16) 

where Ucc is the total cross section for the process e^e~ — >• cc while D^^q^^^J^z) is the respec- 
tive fragmentation function. Instead of the variable ^, its approximate value = p/pmax is 
used below. 

Neglecting a small asymmetry in the angular distribution of Occc baryons, we arrive at 
the following relation between the differential cross section with respect to the transverse 
momentum and the fragmentation function: 

da _ AaccPT } Dc^Q^^^{z)dz ^^^^ 



dpT s z^z^ - Apl/s 



We now compare our numerical results with those obtained according to expression (17) 
with various fragmentation functions. First, we consider the Peterson function [9] 

which is often used in processing experimental data on charmed-hadron production in e'^e~ 
annihilation [14] . Also, this function provides a good approximation to numerical results on 
flgcb baryon production in e+e" annihilation [4]. The best fit to our calculations (dash-dotted 
curve in Fig. 2) corresponds to e = 0.92 . The agreement is clearly poor. The best fits with 
the CoUins-Spiller fragmentation function [10] 

at £ = 3.0 and with the fragmentation function [11] 

D(z) z^'il- z), (20) 

at a = 0.8 are also unsatisfactory. The results of the calculations according to (17) with the 
functions in (19) and (20) are displayed in Fig. 2 (dashed and dotted curves, respectively). 

An acceptable analytic form for our numerical results is provided by the LUND fragmen- 
tation function [12] 

i^(^)^l(l_^)«exp(-^), (21) 

at the parameters a = 2.4 ± 0.2 and c = 0.70 ± 0.03. The corresponding results calculated 
according to (17) are represented by the solid curve in Fig. 2. 

It should be noted that the fragmentation functions (19)-(21), along with the functions 
in (18), were employed by the OPAL Collaboration [15] in processing experimental data on 
S-meson production. 
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1 O'da/dpT, fb/GeV 1 O'dcj/dY, fb 




Fig. 2. Differential cross sections for flccc baryon production in e+e" annihilation at the Z pole 
with respect to the transverse momentum pt (left) and rapidity Y (right). The results of Monte Carlo 
calculations and the errors in them are represented by crosses. The curves in Fig. 2 (left) correspond to 
expression (17) calculated with the fragmentation functions in the form (solid curve) (21) at a = 2.4 and c 
= 0.70, (dash-dotted curve) (18) at £ = 0.92, (dashed curve) (19) at e = 3.0, and (dotted curve) (20) at e 
= 0.8. 
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